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We report on theoretical studies of a recently discovered strong radiation-induced magnetoresis- 
tance spike obtained in ultraclean two-dimensional electron systems at low temperatures. The most 
striking feature of this spike is that it shows up on the second harmonic of the cyclotron resonance 
and with an amplitude that can reach an order of magnitude larger than the radiation-induced re- 
sistance oscillations. We apply the radiation-driven electron orbits model in the ultraclean scenario. 
Accordingly, we calculate the elastic scattering rate (charged impurity) which will define the un- 
expected resonance spike position. We also obtain the inelastic scattering rate (phonon damping), 
that will be responsible of the large spike amplitude. We present a microscopical model to explain 
the dependence of the Landau level width on the magnetic field for ultraclean samples. We find 
that this dependence explains the experimental shift of the resistance oscillations with respect to 
the magnetic field found in this kind of samples. We study also recent results on the influence of 
an in-plane magnetic field on the spike. We are able to reconcile the obtained different experimen- 
tal response of both spike and resistance oscillations versus an increasing in-plane field. The same 
model on the variation of the LL width, allows us to explain such surprising results based in the 
increasing disorder in the sample caused by the in-planed magnetic field. Calculated results are in 
good agreement with experiments. These results would be of special interest in nanophotonics; they 
could lead to the design of novel ultrasensitive microwave detectors. 
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I. INTRODUCTION 

Quantum Hall effect and radiation-matter coupling^ 
are two of the most remarkable topics in Condensed Mat- 
ter Physics. Accordingly, transport excited by radiation 
in a two-dimensional electron system (2DES) has been al- 
ways a central topic in basic and applied research. In the 
last decade it was discovered that when a 2DES in a per- 
pendicular magnetic field {B) is irradiated, mainly with 
microwaves (MW), some striking effects are revealed: 
radiation- induced magnetoresistance {Rxx) oscillations 
and zero resistance states (ZRS)^'-. These remarkable 
effects show up at low B and high mobility samples. It 
is important to achieve a complete understanding of the 
physical mechanisms being responsible of them and not 
only for the basic knowledge purpose but also for the po- 
tential applications in nanoelectronics. Different theories 
have been proposed to explain these effects^"— but the 
physical origin is still being questioned. In the same way, 
a great effort has also been made from the experimental 
sidei^"— . Of course the obtained experimental results 
always mean a real challenge for the existent theoretical 
models. Therefore, a comparison of experiment with the- 
ory could help to identify the importance of the different 
approaches in these theories. Thus, as an example, it 
has been recently published experimental results on the 
dependence of the oscillations with radiation power— i^ 
where a very solid result has been obtained in terms of 
a sublinear relation, similar to a square root. Yet, some 
theories predicted a linear dependence^ between R^x os- 
cillations and radiation power. 

One of the most interesting and challenging experi- 



mental result, recently obtaine d^^'^^ and as intriguing as 
ZRS, consists in a strong resistance spike which shows 
up far off-resonance. It occurs at twice the cyclotron fre- 
quency, w ~ 2w,Mi^, were w is the radiation frequency 
and Wc the cyclotron frequency. The amplitude of such 
a spike is very large reaching an order of magnitude re- 
garding the radiation-induced R^x oscillations. In the 
same experiments it is also reported the behavior of the 
resistance spike with different radiation frequencies, tem- 
perature (T) and radiation power (P). Remarkably, the 
only different feature in these experiments^"''^'^ is the use 
of ultraclean samples with mobility /i ^ 3 x \{f crn^ /V s 
and lower temperatures, T ^ QAK. Yet, for the pre- 
vious, "standard", experiments^!^ the mobility is lower, 
(/i < lO^cmV^s) and T higher, [T > l.OK). Much 
more recently experimental results on this off-resonance 
giant resistance spike have been presented by the same 
authors^. This time they add an in-plane magnetic field 
to the standard experimental set-up. Remarkably, they 
obtain that the radiation-induced resistance oscillations 
and spike present different response to an increasing in- 
plane magnetic field. Resistance oscillations are shifted 
to higher perpendicular magnetic field (B), when the in- 
plane is increased. Yet, the spike approximately keeps 
the same position. Then, the authors suggest that both 
effect are ruled by different physical mechanisms. 

Other strong resistance spikes in radiation-induced 
Rxx oscillations experiments were previously obtained 



around the cyclotron resonance w 
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for higher 



frequencies and power and being explained by resonant 
heating of electrons. Off-resonance phenomena are al- 
ways remarkable phenomena in all branches of Physics. 
They are counterintuitive with respect to the basic reso- 




Energy 



03 CO 
O 



o 



o 

CD 



CD 

E 

CD 
O 
CO 

o. 

CO 




B (a.u.) 



FIG. 1: Upper panel: schematic diagram showing the den- 
sity of Landau states simulated by Lorentzian functions for 
wide and narrow Landau levels, where the width of the states 
is indicated by F. The arrow 1 corresponds to an elastic scat- 
tering process and the arrow 2 to an inelastic one. Lower 
panel: schematic diagram showing the radiation-driven elec- 
tron orbit center displacement (two dimensional electron sys- 
tem) for ultraclean and standard samples. We observe that 
the first one is delayed regarding the latter as being driven by 
radiation of smaller frequency. 



nance process, meaning a challenge to the available the- 
ory. Then, it would be important to fully understand 
these phenomena, not only from the basic standpoint but 
also from the application side. For instance, they could 
lead to the design and development of novel ultrasensi- 
tive photon detectors^ in the microwave and terahertz 
bands, where the technology is not well-developed yet. 

In this article, we theoretically study this radiation- 
induced Rxx spike applying the theory developed 
by the authors, the radiation- driven electron orbits 
mode^'^i^—. According to it, when a Hall bar is il- 
luminated, the electron orbit centers perform a classical 
trajectory consisting in a harmonic displacement along 



the direction of the current at the w frequency. This mo- 
tion is damped by the interaction of electrons with the 
lattice ions and with the consequent emission of acous- 
tic phonons. Thus, the 2DES moves periodically at the 
radiation frequency altering dramatically the scattering 
conditions and giving rise eventually to Rxx oscillations 
and ZRS. We extend this model to a ultraclean sample, 
obtaining that all the scattering conditions are modified. 
Mainly because the LL, which in principle are broad- 
ened by scattering, become very narrow in this kind of 
samples. This implies an increasing number of states at 
the center of the LL sharing a similar energy. In be- 
tween LL, it happens the opposite, the density of states 
decreases dramatically (see upper panel of Fig. 1). There- 
fore, the elastic scattering rate, due to charged impuri- 
ties, increases, meanwhile the inelastic one (phonon emis- 
sion damping) decreases. Thus, the elastic scattering 
time (inverse of elastic scattering rate) decreases reach- 
ing a limiting value of half of the time of a standard 
sample. This specific value is related with the LL degen- 
eracy. For the irradiated electrons, half of the scattering 
time is physically equivalent as being driven by radiation 
of frequency w/2. Accordingly, the cyclotron resonance 
is shifted to a new i?-position around w ~ 'iwc- On 
the other hand, the inelastic scattering decreases and the 
emission of acoustic phonons is less efficient producing 
a bottleneck effect which prevents from releasing the ab- 
sorbed energy to the lattice. Finally the corresponding 
amplitude abruptly increases giving rise to a strong re- 
sistance spike. 

We also present a microscopical theoretical approach to 
the dependence of the LL width on the magnetic field for 
the regime of ultraclean samples. We apply this approach 
in the framework of the general theory of the radiation- 
driven electronic orbits model. We find that this depen- 
dence is very important to explain the experimental shift 
found in the resistance oscillations with respect to the 
magnetic field for this kind of samples^^. We study also 
very recent results on the influence of an in-plane mag- 
netic field on the spike and the corresponding connection 
with the resistance oscillations. We consider that the 
main effect of adding this parallel magnetic field is to in- 
crease the disorder perceived by the electrons in the sam- 
ple. Within our model and among the different sources 
of scattering, inelastic scattering of electrons by inter- 
action with acoustic phonons (emission) is the most di- 
rectly affected. The basic result is a progressive damping 
of the spike and the whole radiation resistance response. 
This is observed in experiments^S. The total scattering 
rate refiects also this increasing disorder which broadens 
the LL; now the LL width depends on both the perpen- 
dicular and in-plane magnetic fields. As feedback effect, 
this will eventually affect the charged impurity scattering 
rate. This feedback and the LL degeneracy cut-off value 
explain the different behavior of the spike and resistance 
oscillations with the in-plane magnetic field. Initially, 
all spike-related phenomena were tried to be explained— 
as unique events with an origin different from radiation- 



induced Rxx oscillations. Yet, in this article we consider X is the center of the orbit for the electron spiral motion: 
that the resistance spike and everything related with it, 
are the outcome of an extreme scenario of the radiation- 
induced oscillations in ultraclean samples and therefore, 
sharing the same physics. 



II. THEORETICAL MODEL 



A. Summary of the radiation-driven electron orbit 
model 
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The radiation- driven electron orbits model, was devel- 
oped to explain the Rxx response of an irradiated 2DEG 
at low magnetic field. We first obtain an exact expression 
of the electronic wave function. Then, the total hamilto- 
nian H can be written as^: 



H 



p2 

2m* 



1 



+ -m*wi{x - Xy - eEdcX + 



1 E'^ 
+ -m*-^ ~ eEo coswt{x - X) - 

— eE() cos wtX 
^ Hi — eEo cos wtX 



(1) 



Eq the intensity for the MW field and Edc is the DC 
electric field in the x direction. Hi is the hamiltonian 
corresponding to a forced harmonic oscillator whose orbit 
is centered at X. Hi can be solved exactl y^^''^- , and using 
this result allows an exact solution for the electronic wave 
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and the main result that we want to point out is that 
^Ar(a;, t) (X 4)n{x — X — Xci{t), t) where (pn is the solution 
for the Schrodinger equation of the unforced quantum 
harmonic oscillator. J„i are Bessel functions and L is 
a classical Lagrangian of an electron in the presence of 
magnetic field. Xd (t) is the classical solution of a forced 
and damped harmonic oscillator: 
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where Eq is the MW electric field. 7 is a damping factor 
for the electronic interaction with the lattice ions giv- 
ing rise to emission of acoustic phonons. Then, the ob- 
tained wave function is the same as the one of the stan- 
dard quantum harmonic oscillator where the center is dis- 
placed by Xci{t). This implies that the electron orbit cen- 
ters oscillate harmonically at w. This radiation — driven 
behavior will affect dramatically the charged impurity 
scattering and eventually the conductivity. 

Next, we calculate the scattering suffered by the elec- 
trons due to charged impurities (elastic) applying time 
dependent first order perturbation theory. Thus, we cal- 
culate the scattering rate^i^i^ between two oscillating 



Landau states (LS), the initial, n, and the final, m: 



1 /■* 

- / < ■^^{x,t)\Vs\^n{x,t) > e'^'dt 



Wn,m = fim — 

a^O dt _ _ 

(5) 

where Vs is the scattering potential for charged 
impuritiei 
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ys = j: 



2Se{q + qa) 



iq- r 



(6) 



S being the surface of the sample, e the GaAs dielec- 
tric constant, and go is the Thomas-Fermi screening 
constant^i^i. After some algebra we arrive at an inter- 
mediate expression for the charged impurities scattering 
rate: 
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r is the LL width and the last part of the expression 
represents the density of final Landau states, pm, in the 



En and Em are the corresponding LL energies for the 
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FIG. 2; a) Calculated irradiated magnetoresistance vs static 
magnetic field for a radiation frequency of / = 103.08GHz, 
different samples (standard and ultraclean) and temperatures. 
For ultraclean sample and a temperature of, T — OAK, we 
observe a intense spike at w « 2wc- Yet, for higher temper- 
ature, T ~ 2.5K, the spike vanishes. For a standard sam- 
ple and low temperature, we obtain a curve with the usual 
radiation-induced resistance oscillations and zero resistance 
states but the spike does not show up. b) Blown-up of upper 
panel for lower magnetoresistance values. We observe that in 
ultraclean samples, radiation-induced oscillations are shifted 
to lower magnetic field positions, compared to standard sam- 
ple case. 



elastic scattering jump. 
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FIG. 3: Calculated irradiated magnetoresistance vs magnetic 
field for radiation frequencies / = 60.0, 103.08, 137.0 and 192 
GHz. For each frequency we obtain a spike at u; « 2wc- The 
spike amplitude increases with the radiation frequency, and 
with the magnetic field as B . This dependence can be shown 
in the inset, where we present the amplitude of the spike vs 
the square of the magnetic field. 



initial and final states respectively. 

Once we know the scattering rate, we consider that 
when an electron undergoes a scattering process jumping 
from the initial state to the final one, it advances an 
average effective distance^ i^^'^^ , AX^"-'^, given by 



AX™'' = AX° + Acosi 



(9) 



where AX'^ is the effective distance advanced when there 
is no radiation field present. Then, following Ridley22, we 
can obtain the expression for the longitudinal conductiv- 
ity (^xx according to the expression. 
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being E the energy. Then, we get to the final expression 
which reads, 
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where ni is the impurities density and Xs 
obtain Rxx we use the relation R^x — —2 
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The radiation power, P, can be related with A through 
the well-known formula that gives radiation intensity / 
(power divided by surface) in terms of the radiation elec- 
tric field Eq: I = ^c£oEq, where c is the speed of light in 
vacuum and ep is the permittivity in vacuum. If we want 
to express only the power in terms of the radiation elec- 
tric field we have to take into account the sample surface. 
In the particular case of GaAs we can readily obtain: 



P — -CCaAs^^oE^S 



(13) 



where ccaAs is the speed of light in GaAs and e is the 
GaAs dielectric constant. Accordingly, 



Ea(xy/P 



(14) 



Then, substituting in the expression of A, we obtain that 
Rxx varies with P following an square root law: 
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Thus, we expect that the Rxx response will grow accord- 
ing to the square root of P, i.e., following a sublinear 
dependence which is finally reflected in the amplitude, 
A(xP°-^. 



B. Far off-resonance position of the giant 
resistance spike 



The precise shape of broadened LL in real systems 
(with disorder) remains, even currently, controversial. 
Usual assumptions are Gaussiai*^"— , Lorentzian^"— 
and semielliptic'^'^, being the two first, the most com- 
monly used. For instance, some published experimental 
results'*^''*-^ indicate that they could be equally described 
or explained assuming either a Lorentzian or Gaus- 
sian profiles. More recently, other works related with 
radiation- induced Rxx oscillations and ZRS, have used 
Lorentzian shapes for LL in the corresponding theoretical 
models^*'*'* In our case, dealing with ultraclean samples 
(very narrow LL), and low B and T, there is no much the- 
oretical or experimental information about the real pro- 
file of the broadened Landau states. We have to consider 
that the samples used in the experiment o^'^'^^ with ex- 
tremely high mobilities have been very recently obtained. 
Then, we are facing a completely brand new scenario re- 
quiring more theoretical and experimental work. There- 
fore, as a first approach and following T. Ihn^, we have 
assumed for the broadened density of Landau states a 
Lorentzian function being the width independent of the 
LL index. 
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If we rewrite eq. (8) considering that £",„ = hwc{m 
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Now, this last expression can be further developed if we 
apply the Poison sum rules^. 
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and then, we can reach for p„ 
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that with En ~ hwdn + 2) ^^ finally obtain, 
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Normally in this calculations if the exponent is not 
very small it is enough to consider the first term s — 1. 
Yet, since we are dealing with high mobility samples and 
very small T the latter approximation can not be applied. 
Then if the total sum is carried out we obtain^: 



that translated into the final density of states, 
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Substituting this result in eq. (7), we arrive at the fi- 
nal expression for the charged impurities scattering rate 
which reads: 
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ing, such as the one caused by charged impurities, most 
processes take place at the center of the LL, sharing ap- 
proximately the same energy (see Fig.l upper panel) and 
then Em — £-«■ This implies that. 
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(25) 



where Wq--^ represents the charged impurities scattering 
rate for a standard mobility sample which is recovered 
for wide F, 
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eJ^ Bra* riiS 



(23) 



As we said above, for ultraclean samples F is very small 
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it turns out 



and for experimental magnetic field; 
that F << hwc- Then, if B increases, the exponent ^^^ 
decreases making Wm,n continuously to increase. How- 
ever, there exists a cut-off defined by the LL degeneracy: 



LLdegen. 



2eB 



(24) 



In this regime of narrow F, in an interval energy of hwc 
we must have, in average, a number of states given by 
the LL degeneracy and therefore: 



2eB 1 ^ 2eB 



(26) 



This is only fulfilled if. 
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In a regime of narrow LL the density of states pm can 
be approximated considering that in an elastic scatter- 



which is a remarkable result because translated into the 
expression of eq. (22) we obtain: 
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Then, the charged impurities scattering rate increases 
with B till a cut-off value of twice the corresponding of 
a standard sample. In terms of charged impurities scat- 



tering time. 



^uclean — ^ '^ 



(29) 



where r = I/Wq is the time of a standard sample and 
Tuciean corresponds to a ultraclean one. Accordingly, the 
ultraclean scattering time, Tudean = pfr — turns out to be 
twice shorter than the standard and then, the scattering 
event is twice faster. 

Within the framework of our theory'* , this implies that 
during the time Tuciean and compared to the standard for 
the same B, the ultraclean 2DES appears to be displaced 
by radiation a smaller distance. In other words, in terms 
of scattering, the ultraclean harmonic motion (electron 
orbit center displacement) is perceived as delayed regard- 
ing the standard, as if electrons were driven by radiation 
of smaller frequency (see Fig. 1 lower panel) , more pre- 
cisely, of half frequency. Then, the radiation electric field, 
Ew, is perceived as. 



E„ 



En COS —t 

2 



(30) 



by the ultraclean 2DES. The conclusion is that in ultr- 
aclean samples, during scattering, electrons "feel" radi- 
ation with half of the real frequency. Applying next the 
theory^, we reach an expression for R^x in the ultraclean 
case: 
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with the lattice ions, damping the electronic orbits mo- 
tion and releasing radiation energy in form of acoustic 
phonons. Therefore, 7 is given by^: 
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(32) 
where Wac is the frequency of acoustic phonons for exper- 
imental parameters^i^ and the last term represents the 

density of final Landau states. Following similar steps as 
before we obtain the expression: 
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where S is the acoustic deformation potential, p the mass 
density, fcs the Boltzman constant, T the temperature, 
vs the sound velocity and < z > is the effective layer 
thickness. When the total sum inside brackets is carried 
out we obtain^. 



According to it, now the resonance in R^x will take place 
at w sa 2wc, as experimentally obtained^iiS^. 



C. Intensity of the resistance spike 



E 



exp 



hwc 



cos 



2'KshWa 
hWr 



The intensity of the Rxx spike will depend on the rel- 
ative value of the frequency term, (w^ — (y)^)? &n.d the 
damping parameter 7 in the denominator of the latter 
Rxx expression. When 7 leads the denominator the spike 
is smeared out. Yet, in situations where 7 is smaller than 
the frequency term, the resonance effect will be more vis- 
ible and the spike will show up. As we explained above, 
the parameter 7 represents the interaction of electrons 
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Then, substituting in 7 we finally obtain. 
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In order to have an approximate but more compact 
version of eq. (27) we can consider average values for 
the cyclotron energy and for the acoustic phonon energy 
in the case of experimental parameters. Thus, for the 
magnetic fields swept in the experiments we can take in 



average that the cyclotron energy is hwc 
and the phonon acoustic energy is, hwac 



2xlO-*meV, 
1 X 10"* meV. 
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Then, for those values the cosine term, 
2nhwac 
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(36) 



that gives an expression for the damping parameter 
which reads 
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According to this last expression, for increasing B the 
term inside brackets decreases and as a result, the pa- 
rameter 7 will get smaller, making increasingly difficult 
the damping by acoustic phonon emission and the re- 
lease of the absorbed energy to the lattice. Therefore, 
we have a bottleneck effect for the emission of acous- 
tic phonons. Now it is possible to reach a situation 
where {w1 — (y)^)^ ^ 1^ making visible a resonance 
effect and, therefore, giving rise to a strong resonance 
peak at w sa 2wc- For GaAs and standard experimen- 
tal parametersSsi, we obtain that 7 ~ 7 — 10 x lO^^s"^ 
and that in a ultra-clean regime^^i^ it decreases till 
~ 3.5 X 10^^ . This fixes a lower cut-off value for the 
radiation frequency where the resistance spike could be 
observed. According to our calculations this would be 
around / = w/27r sa 40 — 4:5GHz. Experiments have 
also pointed out that at lower frequencies is much more 
difficult to observe resistance spikes. 



D. Variation of F with the magnetic field for 
ultraclean samples 

One important issue which plays a key role in the 
model presented above is the dependence of T on the 
magnetic field. Previous experimental and theoretical 
works^^"— report that for lower B (but higher than 1 
T), r varies with the magnetic field as ~ y/B. Yet, for 
increasing B, F becomes independent. For our brand 
new scenario of extremely high mobility samples at very 
low B, {^ 0.1 T), and very low temperature, (~ 0.4 K), 
there is no much information, experimental or theoreti- 
cal, about the actual dependence of F on B. Therefore, 
and as a first approach, we have developed a simple mi- 
croscopical description, based in scattering parameters, 
of this dependence which starts from previous analytical 
expressions 



33.35.37.42.45.48 



Thus, following T. Ihnsa, we 



begin with the expression which relates F and B, 
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where tq is the zero magnetic field quantum scattering 
lifetime. 

For the experimental spikea^ii^, the magnetic field 
used is very low, around B '^ O.IT. Then, for this regime 



we can assume that the variation of F with B is closer to 
a straight line than to a square root and proceed to lin- 
earize the latter expression. Thus, we expand the square 
root by a Taylor series around B — 0.1 T, till the first 
derivative term. Then, we can write. 



T{B) = r(o.i) + {B- o.i)r'(o.i) 



(39) 



and finally obtain the expression which gives the depen- 
dence of F on _B for very low B 



r{B) oc 




(40) 



where p — > 1. p is a phenomenologically introduced pa- 
rameter which takes into account that the dependence 
of F on _B does not follow strictly a straight line but it 
is close to it. In our case we have assigned p ~ 0.90. 
This expression implies that, for extremely high mobil- 
ity samples and very low B, the Landau states broaden 
almost linearly with the magnetic field affecting all scat- 
tering processes such as elastic impurity scattering and 
inelastic phonon scattering. This dependence of F on i? 
has been used throughout the above theoretical model on 
spikes. 



E. Influence of an in-plane magnetic field 

The effect of an in- plane magnetic field, (-B[|), on 
the radiation-induced resistance oscillations was al- 
ready studied and published^? , and experimental results 
showed that the main effect was a progressive damping 
of the whole resistance response as i3|| increased. Subse- 
quent theoretical results^'', confirmed and explained the 
surprising damping in the framework of the radiation- 
driven electrons orbits model. Accordingly and as we said 
above, in an irradiated 2DES in the presence of a perpen- 
dicular B, electronic orbits are forced to move back and 
forth, oscillating harmonically at the frequency of radia- 
tion and with an amplitude proportional to the radiation 
electric field. In their radiation-driven motion, electrons 
interact with the lattice ions being damped and produc- 
ing acoustic phonons. The presence of B|| imposes an 
extra harmonically oscillating motion in the z-direction 
enlarging the electrons trajectory in their orbits. This 
would increase the interactions with the lattice making 
the damping process more intense and reducing the am- 
plitude of the orbits oscillations. Therefore, the effect 
of the presence of B|| is to increase the disorder in the 
sample from the electrons perspective. 

The Hanfiltonian for electrons confined in a 2D system 
(x-y plane) by a potential V{z) and subjected to a to- 
tal magnetic field Bt = {Bx,0,B), (i?|| = B^) is given 
considering the previous hamiltonian of eq. (1) by: 



H = Hi- eEo cos wtX + H^ 



(41) 
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X 
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FIG. 4: Calculated irradiated magnetoresistance vs static 
magnetic field for different radiation power (P) and a fre- 
quency / = WSmGHz. P varies from 3.8 mW to 0.1 mW, 
including P = 2.7, 1.7, 1.0 and 0.4 mW. For decreasing P, 
the radiation-induced magnetoresistance response decreases 
as well, and for P — O.lmW we obtain a response close 
to darkness. In the inset we present the spike amplitude 
vs P. We fit the data obtaining a sublinear P-dependence, 
Rxx C)c P" where a is close to 0.5. 



If we consider a parabolic potential for Vz 



Viz) 



-m WnZ 



(42) 



Hamiltonian 


H 


J can 


be written as 




Hz 


= 


Pz 


, * / 2 , 

+ 2™ (^x + 


wl)z 


2m* 




= 


P! 


+ ^m*n^z^ 






2m* 





We have used the the Landau gauge for B^ 
{0,—zBj;,0), and 



Wx 



eBx 

m* 



(43) 



Ab. - 



(44) 



The use of a parabolic potential for Vz allows that the 
corresponding time-dependent Schrodinger equation of 
H can be readily solved. We obtain the wave functions 
of two harmonic oscillators, one in the x direction and 
the other in the z-direction; 



In a semiclassical approach the electron is subjected si- 
multaneously to two independent harmonic motions with 
a trajectory depicted in Fig. 9: As we have indicated 
above, the presence of in-plane B alters the electron tra- 
jectory in its orbit increasing the frequency and the num- 
ber of oscillations in the z-direction. Now the frequency 
of the z-oscillating motion is Vt > wq. This makes longer 
the electron trajectory increasing the total orbit length 
and eventually the damping. This increase in the orbit 
length is proportionally equivalent to the increase in the 
number of oscillations in the z-direction. Thus, we intro- 
duce the ratio of frequencies after and before connecting 
Bx as a correction factor for the damping factor 7. The 
final damping parameter 7/ is: 



n 



7/ = 7X — = 7X 

Wo 




7X '1/1 + 



h 



2\ 2 



(46) 



where zq is the effective length of the electron wave func- 
tion when we consider a parabolic potential for the z- 
confinemen^i^, 



zo 



mrwo 



(47) 



The expression of 7/ shows that one of the main results 
of the presence of Bx is an increase in the damping pa- 
rameter due to the interaction with lattice ions. There- 
fore, the whole resistance response to radiation, (spike 
and oscillations) will be increasingly damped. The sec- 
ond important effect is reflected in the total quantum 
scattering rate I/tq. According to the Matthiessen rule 
the total scattering rate can be expressed as the sum of 
the different individual scattering sources, — — V, — 
and obviously one of them is the acoustic phonon scat- 
tering rate. The increase in the phonon scattering rate 
(7 = 1/tqc) due to Bx will eventually affect the total 
scattering rate that will increase as well, making it Bx- 
dependent, ^ tg \ - Considering the Matthiessen rule and 
the ratio of frequencies — affecting 7, we can readily ob- 
tain the expression, 

1 



ro{B:, 



1 

- + 7 
7-0 


"-1 

Wq 






1 

- + 7 

To 


[h 




-1 



'^Tix.Z.t) (X (j)N [{x~ X -Xcl{t)),t]4){z) 



(45) 



(48) 
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And now eq. (40) can be rewritten taking into account 
the last expression, 



T{B,B^) 



1 



n?e 



2yoTV 27rm*ro(5^ 



-BP (49) 



Now r depends on both magnetic fields the perpendicular 
B and the in-plane B^, 



III. RESULTS 

In Fig.l in the upper panel, we present a schematic 
diagram showing the density of Landau states for wider 
and narrow F, (standard and ultraclean samples respec- 
tively). Therefore, we can observe the two different 
regimes; for a standard sample where hwc < F and for a 
ultraclean one where hwc > F. The broadened Landau 
states have been simulated by Lorentzian functions. The 
arrow 1 corresponds to an elastic scattering process (re- 
mote charged impurity scattering) and the arrow 2 to an 
inelastic one (acoustic phonon scattering). In the lower 
panel of Fig. 1, we present the radiation-driven electron 
orbit displacement for ultraclean and standard samples. 
We observe that the first one is delayed regarding the 
latter as being driven by radiation of smaller frequency. 

In Fig. 2a, we present calculated irradiated R^x vs 
static magnetic field for a radiation frequency oi f — 
103.08GHz. The presented results correspond to stan- 
dard and ultraclean samples and high and low temper- 
atures. For the ultraclean sample and T = O.AK, we 
obtain a strong spike at w « 2wc as in experiment o^*^'^^ . 
Increasing temperature for the same sample, (T ~ 2.5K) 
the spike vanishes but the radiation-induced oscillations 
still remain but with lower intensity, as expected. Finally, 
for a standard sample, we obtain the usual radiation- 
induced Rxx oscillations and ZRS. In Fig. 2b we present 
a blown up of the lower values of Rxx of the upper panel. 
Thus, we can contrast the curve of the ultraclean sample 
with the standard, observing a shift of the Rxx oscilla- 
tions to lower magnetic fields. 

In our calculations we have used a quantum lifetime at 
zero magnetic field tq ~ 2 x 10~^^ g ^24. 25,51 f^j. ^^^ ^■^^^._ 
aclean Rxx curve, which gives us, according to eq. (32), 
a LL width of F - 5. x lO-'^B^ eV. For the standard 
curve, we have used a shorter quantum lifetime,— i^^— 
To - 10^12 s^ giving a broader F, F ~ 25.10~'^BP eV. 
The quantum life time defines the quality of the sample; 
the longer this time, the higher the quality and narrower 
F. Smaller F, according to eq. (22), will give us a longer 
impurity scattering time. Then, the perceived radiation 
frequency will be smaller, shifting the Rxx oscillations 
to lower B as observed in ultraclean samples. On the 
other hand, with a shorter quantum life time (low qual- 
ity samples or standard), F is large and the term inside 
brackets in eq. (22) tends to 1. Thus, we recover the 
usual positions for radiation-driven Rxx oscillationsii. 
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FIG. 5: Calculated irradiated magnetoresistance vs magne- 
toresistance for different temperatures and a frequency / = 
lOS.OSGHz. Temperature varies from 0.4 k to 4.0 K. We ob- 
serve a clear decrease of the spike and Rxx oscillations for 
increasing T. The T-dependence is explained with the damp- 
ing parameter 7 which represents the interaction of electrons 
with acoustic phonons. In the inset we present the ampli- 
tude of the spike vs T. We fit the data obtaining a relation 
Rxx oc T~^ (hyperbole), as expected from the model. 



In Fig. 3 we present calculated irradiated Rxx vs mag- 
netic field for different radiation frequencies / = 60.0, 
103.08, 137.0 and 192 GHz. For all cases we obtain a 
clear spike at w « 2wc. We observe that the spike am- 
plitude increases with the frequency and as a result with 
the magnetic field. This is because of the presence of 
the frequency term, (w^ — ^ )^ in the denominator of 
the Rxx expression, and then, the higher w, the higher 
the magnetic field where the resonance takes place. In 
the inset of the figure we observe that the amplitude of 
the corresponding spike grows with the magnetic field 
following a B^ law. This expected from our model^ be- 
cause according to it, the radiation-dependent part of 
the Rxx depends on the magnetic field as B^ . Therefore, 
we predict for the amplitude of the "w « 2wc spike", 
a square law dependence on B"^ that would need to be 
experimentally confirmed. The experiments that discov- 
ered the Rxx spike showed that the amplitude increases 
with the frequency. Yet, they do not present any law or 
fit showing the increase rate of the spike amplitude with 
B. 

In Fig. 4 we present calculated irradiated Rxx vs mag- 
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FIG. 7: Calculated values for A versus B. For B ~ the term 
A — !> 1 because the absolute value of the exponents is very 
large. Then, as B increases the exponents decrease and A 
increases till the cut-off value of 2. This makes the impurity 
scattering rate continuously to increase with B from A ~ 1 
to A ~ 2, or impurity scattering time to decrease. This effect 
is perceived by the electrons as if radiation had a decreasing 
frequency. Eventually, the back and forth motion of electrons 
in their orbits is performed at a decreasing frequency too. 



FIG. 6: Normalized B~ plots of calculated data on radiation 
induced resistance oscillations of Fig. 2 for standard sample 
(6a) and for ultraclean (6b). S's are the oscillatory periods in 
B~^ for each case. 



netic field for different radiation power (P) and frequency 
/ = 103.08GHz . P varies from 3.8 mW to 0.1 mW, in- 
cluding P = 2.7, 1.7, 1.0 and 0.4 mW. Decreasing P, 
the radiation-induced Rxx response, (spike and oscilla- 
tions), decreases as well, and for P = O.lmM^ we obtain 
nearly the darkness result. In the inset we present the 
amplitude of the R^x spike vs P. We fit the data ob- 
taining a sublinear P-dependence, Rxx oc P" where a is 
close to 0.5, as expected, and explained in terms of our 
model: Eq cx vP => Rxx oc vP and in agreement with 
experimental^! and theoretical^ results. 

In Fig. 5 we present calculated irradiated magnetore- 
sistance vs magnetic field and frequency / = 103.08GHz, 
for differente temperatures (T). T varies from 0.4 to 4.0 
K. We observe a clear decrease of the spike for increasing 
T, till it is smeared out. The T-dependence, according to 
the model, is explained with the damping parameter 7. 7 
is linear with T'^i^^^ and Rxx oc 7~^. Thus, an increasing 
T means an increasing 7 and smaller both spike and Rxx 




5 6 7 8 9 10 
8 



FIG. 8: Calculated irradiated Rxx vs e for a frequency of 
/ — 103.08 GHz. We observe a qualitatively similar shift as 
in experimen1i24. Yet, quantitatively speaking the calculated 
shift is larger than the obtained in the experimenli^. 
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♦ B, 
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(c) 
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♦ B. 
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FIG. 9: Schematic diagrams showing the semiclassical de- 
scription of electron trajectories in 2D systems. In the three 
of them there is present a perpendicular magnetic field Bz 
that explains the circular trajectories, a) 2DES {x — y plane) 
and Bz- b) 2DES, B^, a parabolic potential T4 and an in- 
plane magnetic field i3||. c) Same as b) but the intensity of 
-B|| is bigger. 



oscillations. When the damping is strong enough (higher 
T) all Rxx response collapses giving a final result close 
to darkness. The inset shows the relation R^x oc T~^ 
(hyperbole), as expected from the model. 

In Fig. 6 we present normalized B~^ plots of calcu- 
lated data of radiation induced Rxx oscillations of Fig. 2, 
(/ = 103.08Gi?z), for standard (6a) and ultraclean (6b) 
samples. (5's are the oscillatory periods in B~^ for each 
case. In the standard case we observe a nearly perfect 
periodic curve with respect to B~^ , as experimentally 
observed^. This is explained with the linear dependence 
on B of the charged impurity scattering rate (see eqs. 
(22) and (23)) when the brackets term tends to 1 (stan- 
dard samples). This implies that the charged impurity 
scattering time depends inversely on B, giving the peri- 
odic curve of Fig. 6a {Rxx oc A cos wt) . In the ultraclean 
case, we observe an oscillatory behavior of Rxx versus 
B^^ but it is not perfectly periodic. Instead, we observe 
a shift of the Rxx oscillations regarding the vertical lines 
defined by the period. Peaks and valleys are increasingly 
shifted. As B increases, the distance between subsequent 
peaks (or valleys) increases. 

The origin of this effect is in the exponents of eq. (22) 
which depend on B. The dependence on B of the charged 
impurity scattering rate comes not only from the linear 
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FIG. 10: Calculated irradiated Rxx versus B for a radiation 
frequency oi f — 104.88 GHz and for different values of the 
in-plane magnetic field Bx- Bx goes from OT till 0.8T in 
intervals of O.IT. We first observe a progressive damping of 
the resistance spike and resistance oscillations. For approxi- 
mately Bx ~ 0.6 T the spike has been completely removed. 
We observe as well that starting from Bx ~ 0.5 and for larger 
Bx the main Rxx associated with the spike shifts to higher B. 
However, we observe that from Bx = T till Bx = 5 T, the 
Rxx spike remains approximately in its position, although its 
intensity is progressively smaller. The green arrows highlight 
this effect. 



term of Wq of that equation, but also from the exponen- 
tials in the brackets term what we call A, 



1 + exp 


hivc 


1 — exp 


hwc 



(50) 



The variation of A with B is presented in Fig. 7. For 
B close to 0, A — >■ 1 because the absolute value of the 
exponents is very large. Then, as B increases the ex- 
ponents decrease and A continuously increases till the 
cut-off value of 2. This makes the charged impurity scat- 
tering rate continuously to increase with B from A ~ 1 
to A ~ 2, (or impurity scattering time to decrease). This 
effect is perceived by the electrons as if radiation had a 
decreasing frequency. Eventually all of this is reflected in 
the obtained Rxx where the distance between Rxx peaks 
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FIG. 11: Calculated irradiated Rxx vs e for a frequency of 
/ = 103.08 GHz. In a) e goes from 1.5 to 3.5. As in Fig. 10 
we observe the spike keeping the same position till Bx = 5 
T. Starting from Bx = 6, the main resistance peak associated 
with the spike shifts to higher B as Bx increases. In b) e goes 
from 2 to 7. We observe the shift, also to higher B, of the 
other resistance peaks. The shift is represented by A^,. Yet, 
in apparent contradiction, the corresponding spike stands still 
for the same Bx values as observed in a). 



(or valleys) gets larger. This is what explains the ob- 
tained shift in Fig. 6b. 

From the experimental standpoint a similar shift can 
be found in Fig. lb of reference [24] . In this Fig. it is 
represented the difference between resistivity under radi- 
ation minus dark resistivity versus the parameter e = -^ 
and obtained an oscillation shift of 0.25 of the periocf. 
According to the experimental results^!, the shift is in- 
creasingly larger for decreasing B. To contrast our cal- 
culated results on this shift with experiment we present 
in Fig. 8 calculated Rxx vs e. We obtain a qualitative 
agreement because we obtain a similar variation of the 
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FIG. 12: Calculated values of A term versus B for same val- 
ues of Bx as in Fig. 10, and Bx ~ 1.1 T. The interpretation 
is similar to Fig. 7. The presence of an increasing Bx, gives 
rise to larger F which as a result makes A to decrease. Then, 
the impurity scattering time decreases continuously being re- 
flected in the perceived radiation frequency which decreases 
too. The final outcome is that Rxx oscillations shift to higher 
B as Bx increases. This is a general behavior for all values of 
Bx- We can see in the figure that from Bx = T till Bx = 5 
the threshold of A ~ 2 is achieved in the sweeping B before 
the resonance condition of u; = 2wc or around it. Then, when 
reaching this condition, the resonance effect can take place 
and an intense Rxx spike will rise in the same position for all 
of them, 



oscillations shift with B. Yet, quantitatively speaking we 
obtain a larger value, reaching 0.5 of the period. The rea- 
son for this quantitative discrepancy could be explained 
by the simple microscopical model used to describe the 
dependence of F on _B for ultraclean samples. On the 
other hand, the model is able to explain the existence of 
the shift, its variation with B and the connection with 
the resistance spike. 

In Fig. 9, we present schematic diagrams showing the 
semiclassical description of electron trajectories in 2DES 
(x-y plane), under a perpendicular magnetic field in the 
z-direction {Bz = B). We present three cases. In a) 
we have the basic situation of only a 2DES and B^- In 
the case b) we add to a) a parabolic potential Vz and 
an in-plane magnetic field Bn . And finally in the case c) 
we have the same as b) but the intensity of i3|| is larger. 
In a semiclassical approach the electron dynamics under 
simultaneously these two magnetic fields can be inter- 
preted as being subjected to two independent harmonic 
motions with trajectories depicted in Fig. 9. When one 
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considers only Bz , the electron performs a circular move- 
ment in the x — y plane, (see Fig. 9a). In Fig. 9b, we add 
a parabolic potential in the z direction and we introduce 
Bu , then the electron trajectory is circular in the plane 
and at the same time is oscillating in z. In Fig. 9c Bu 
increases and the number of oscillations in z direction 
increases too. 

In Fig. 10 we present calculated irradiated Rxx ver- 
sus B for a radiation frequency of / = 104.88 GHz and 
for different values of the in-plane magnetic field Bx ■ Bx 
goes from T till 0.8 T in intervals of 0.1 T. We first 
observe a progressive damping of the Rxx spike and os- 
cillations. For approximately Bx = 0.6 T the spike has 
been completely removed. We observe as well that for 
Bx > 0.5, the main Rxx peak associated with the spike, 
shifts to higher B. As we have previously explained, the 
main effect of Bx is an increase of disorder in the sample. 
In our approach this is reflected in a stronger interaction 
of electrons with the lattices ions giving rise to a more in- 
tense emission of acoustic phonons. Then, the absorbed 
radiation energy can be more efficiently released and the 
resistance response tends to progressively collapse. This 
is similar to the effect of an increasing lattice tempera- 
ture as previously presented (see Fig. 5). Remarkably 
we observe that from Bx = T^ till S^; = 5 T, the Rxx 
spike remains approximately in its position, (see green 
arrows in Fig. 10), although its intensity is progressively 
smaller. 

In Fig. 11a, we present similar information as in Fig. 
10 but this time versus e. In Fig. lib, we present a blown 
up of 11a for e values from 2 to 7. We observe that mean- 
while the spike does not move when we sweep B^ from 
Bx — Q T! till Ba; = 5 T, the corresponding resistance 
oscillations present a clear shift to higher B. One could 
think in principle that such different Bx responses could 
indicate different physical origin. However following our 
theory, they are totally connected. And the explanation 
has to do, as before, with the variation of F with the 
magnetic field. The explanation can be obtained from 
Fig.l2. 

In Fig. 12 we present the A term versus B for same 
values of Bx as in Fig. 10, and adding Bx = 1.1 T. In 
principle, the interpretation is similar to Fig. 7. Now 
the presence of an increasing B^,, gives rise to larger 
F which as a result makes A to decrease. Then, the 
impurity scattering time decreases continuously being 
reflected in the perceived radiation frequency which 
decreases too. The final outcome is that Rxx oscillations 
shift to higher B as Bx increases. The effect is similar 
as having increasingly dirtier samples. This is a general 
behavior for all values of Bx ■ The only difference is that 
for increasing B^,, A grows slower versus B. Yet, we 
can see in the figure that from B^, = T till Bx — b 
the threshold of A ~ 2 is achieved in the sweeping B 
before the resonance condition of w; = 2wc or around 
it. Then, when reaching this condition, the resonance 
effect can take place and an intense Rxx spike will rise 
in the same position for all of them, as experimentally 



obtained^. Therefore, we obtain some Bx values, 
[Bx < 5) where constant spike position coexist with a 
shift of Rxx oscillations to larger B. For the Bx values 
that reach A ~ 2 after the resonance condition, the 
spike apparently would rise at higher B but the intense 
damping effect done by Bx makes it to vanish and can 
not be observed any longer. However, the shift of Rxx 
oscillations still can be observed and an example is the 
shift of the main Rxx peak associated with the spike. 



IV. CONCLUSIONS 

In summary, we have theoretically studied the recently 
discovered intense radiation-induced magnetoresistance 
peak obtained in ultraclean 2DES. The most remarkable 
feature of such a peak is that it shows up at w ~ 2wc and 
with an amplitude of an order of magnitude larger than 
the standard radiation-induced magnetoresistance oscil- 
lations. We apply the radiation-driven electron orbits 
model and we calculate the charge impurity elastic scat- 
tering rate, which increases in ultraclean samples, to ob- 
tain the novel resonance peak position (far off resonance). 
We obtain the inelastic scattering rate by phonon damp- 
ing, which decreases in ultraclean samples, showing that 
it is responsible of the large peak amplitude. We present 
a microscopical model to explain the dependence of the 
LL width (r) on the magnetic field for ultraclean sam- 
ples. We find that this i?-dependent variation of F is 
essential to explain the experimental shift found in the 
resistance oscillations. Accordingly, the shift and the re- 
sistance spike are physically connected. 

We have studied also very recent results on the influ- 
ence of an in-plane magnetic field on the spike. Accord- 
ing to them radiation-driven resistance oscillations and 
spike offer different behavior when this magnetic field is 
increased, as if they had different physical origin. The 
same model allows us to explain such surprising results. 
We conclude that the role of the in-plane magnetic field 
is mainly to increase the disorder making the ultraclean 
sample to behave progressively as a dirtier sample and 
affecting the total scattering rate and the width of the 
LL. Now, the simultaneous dependence of this width on 
B and Bx explains the apparent contradiction. Another 
important physical effect observed in experiments^'* for 
the first time in ultraclean samples is an important de- 
crease in the measured resistance for very low magnetic 
fields. This effect is obtained without radiation and is 
described in the experiments as a pronounced negative 
magnetoresistance. Yet, as the magnetic field increases, 
the usual values of standard magnetoresistance are re- 
covered. This effect, closely related with the quality of 
the sample, is the main topic of a future work. Calcu- 
lated results are in good agreement with experiments. 
These results are remarkable also from a technological 
standpoint. For instance in nanophotonics where they 
could lead to the design and development of ultrasen- 
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sitive photon detectors in the microwave and terahertz 
bands where the technology is not mature yet. 
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